Master Operators Govern Multifractality in Percolation 
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Using renormalization group methods we study multifractality in percolation at the instance of 
' noisy random resistor networks. We introduce the concept of master operators. The multifractal 

moments of the current distribution (which are proportional to the noise cumulants (x, x') of 
the resistance between two sites x and x' located on the same cluster) are related to such master 
' operators. The scaling behavior of the multifractal moments is governed exclusively by the master 

operators, even though a myriad of servant operators is involved in the renormalization procedure. 
We calculate the family of multifractal exponents {tpi} for the scaling behavior of the noise cu- 
mulants, c£ (x,x') ~ \x — x'\^ 1 , where v is the correlation length exponent for percolation, to 
QJ • two-loop order. 
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PACS numbers: 64.60.Ak, 05.40.-a, 72.70.+m 

Percolation |l| is a leading paradigm for disorder. Though it represents the simplest model of a disordered system, 
^ h has many applications, e.g., polymerization, porous and amorphous materials, thin films, spreading of epidemics 
etc. In particular the transport properties of percolation clusters have gained a vast amount of interest over the 
last decades. Random resistor networks (RRN) are the most prominent model for transport on percolation clusters. 
Here we discuss RRN in the context of multifractality We propose the concept of master operators ('operator' 
in the sense of 'composite field'). Each moment of the multifractal measure (the square of the bond currents) has a 
master operator as field theoretic counterpart. These master operators are highly and dangerously irrelevant in the 
renormalization group sense. Therefore, each master operator needs in general a myriad of other irrelevant operators 
1 1' for renormalization. However, the renormalization of these servant operators does not induce their master. It follows, 
that the servant operators can be neglected in determining the scaling index of their master operator, i.e., one is 
' spared the computation and diagonalization of gigantic renormalization matrices. 

Systems where multifractality has been observed include, beside RRN, turbulence, diffusion near fractals, electrons 
00 ■ in disordered media, polymers in disordered media, random ferromagnets, and chaotic dissipative systems. Multifrac- 
(T) , tality in RRN is exhibited by the moments of the current distribution measured between two connected terminals x 
and x' . The multifractal moments of the current distribution are related to the noise cumulants (x,x') of the 
resistance by Cohn's theorem jj). At criticality the noise cumulants scale as {x,x') ~ \x — x'\~ x ", where the x n 
constitute an infinite set of exponents which are not related to each other in a linear way as commonly occurs in 
■ critical phenomena under the name of gap scaling. The multifractality manifests itself in the nonlinear dependence 
of x n on n. The existence of the set {x n } was proposed by Rammal et al p|. A set of exponents {C271} equivalent to 
{— x n v}, where v is the correlation length exponent for percolation, was also proposed by Arcangelis et al. [0. 

The field theoretic description of multifractality in RRN was pioneered by Park, Harris and Lubensky (PHL) ||. 
Based on an approach by Stephen [Q they formulated a DxE-io\d replicated Hamiltonian to study the effects of noise in 
RRN. The contributions to the Hamiltonian leading to multifractal behavior contain powers of replica space gradients 
O ■ analog to powers of real space gradients, which were accounted for as an origin of multifractality by Duplantier and 
Ludwig ||. PHL introduced a set of exponents {ip n } identical to the set {— x n i>} and calculated it to first order in 
e = 6 — d, where d denotes the spatial dimension. 

Consider a d-dimensional lattice, where bonds between nearest neighboring sites i and j are randomly occupied 
with probability p or empty with probability 1 — p. Each occupied bond b = (i, j) has a conductance ov Unoccupied 
bonds have conductance zero. The bonds obey Ohm's law cry (Vj — Vi) = where Zy is the current flowing 
through the bond from j to i and Vi is the potential at site i. The <ib are equally and independently distributed 
random variables with distribution function /, mean er, and higher cumulants A^™- 2 ) satisfying A^") <C ~5 n ■ The 
noise average is denoted by {• • = J Y\ b da^f (<Tfc) • • • and its nth cumulant by {• • - n }^\ Both kinds of disorder, 
the random dilution of the lattice and the noise, influence the statistical properties of the resistance R(x, x') of the 
backbone between two sites x and x' . They are reflected by the noise cumulants 
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{x,x') = (x{x,x>){R{x,xT}f) c /{x{x,x')) c , (1) 
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where (■ • -)c denotes the average over the configurations C of the randomly occupied bonds and x[x, x ') is an indicator 
function which is is unity if x and x' are connected and zero otherwise. is related to the 2nth multifractal moment 
of the current distribution via Cohn's theorem [| (cf. @), 



C&° (x,x')=v n {x(x,x , )J2 (h/lf l I (X(x, x')) 
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where v n — {(Spi,) 11 }^ is the nth cumulant of the deviation dp^ = pi, — P of the bond resistance pi, = a h 1 from its 
average p. 

Our aim is to determine Hence, the task is to solve the set of Kirchhoff's equations and to perform the 

averages over the diluted lattice configurations and the noise. It can be achieved by employing the replica technique 

0. PHL introduced DxE-fold replicated voltages V x -^V x— [Vx™"^) and ip~(x) = exp ( i A ■ Vx), where 

V / a,j3— 1 A \ / 

A • V x = Efii X^vt^ and aVo. The corresponding correlation functions are defined as 



G ( x, x': A I = lfm 
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is the power dissipated on the backbone and Z is the usual 



where P ({v}) = £££ 1 E W ><#(^ 
normalization. 

The integrations in Eq. (|^) can be carried out by employing the saddle point method, which is exact in this case. 
The maximum of the integrand is determined by the solution of Kirchhoff's equations and thus 



g(x,x';\) = {ti | exp 
On defining K x (a) = Ef=i [E£=i (A^) 2 ]' one obtai 



i? (/3) {x,x') 
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G ( x f x'] A ) = ( exp 
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i.e., G represents a generating function for C R 

To guarantee that the limit lim^^o Z DE is well defined one switches to voltage variables 9 and current variables 
A taking discrete values on a .DX-E-dimensional torus. Upon Fourier transformation in replica space one introduces 
the Potts spins |§ (x) = (2N'iy DE J2~ « exp ( i A • 9 ) ^«(a;) = « - (2My DE subject to the condition 

9 \j^0 \ /A 9 , 9 x 

^„ (x) = 0. The effective Hamiltonian 
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can be expanded in terms of the Potts spins as 
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The Fourier transform K(\^oiK(0^m turn has the expansion 
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with t, it) ~ ct -1 , and ~ /a 21 being expansion coefficients. We display only the leading parts of the expansion. 
The terms proportional to KA \ J are retained because they are required in extracting the multifractal moments 
from G(x, x'; A ). 

We proceed with the usual coarse graining step and replace the Potts spins (x) by order parameter fields y(x, ) 

which inherit the constraint ^<-><^(x, J = 0. We model the corresponding field theoretic Hamiltonian TL in the 

spirit of Landau as a mesoscopic free energy from local monomials of the order parameter field and its gradients in 
real and replica space. The gradient expansion is justified since the interaction is short ranged in both spaces for large 
a. Purely local terms in replica space have to respect the full S^ 2 m) de Potts symmetry. After these remarks we write 
down the Landau-Ginzburg- Wilson type Hamiltonian 
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In Eq. ( pL0|) we have neglected terms which are irrelevant in the renormalization group sense, r, w and m are now 
coarse grained analogues of the original coefficients appearing in Eq. (||). Note that H reduces to the usual (2M) DE 
states Potts model Hamiltonian by setting vi = and w = as one retrieves purely geometrical percolation in the 
limit of vanishing vi and w. 

Now we analyze the relevance of the vi . A straightforward scaling analysis reveals that 



Cf ((x,x');t,™,M) 
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where // is a scaling function. Note that the coupling constants Vk appear only as Vk/w k . Dimensional analysis of 

the Hamiltonian shows that w A ~ l^ 2 and VkKk( A ) ~ /i 2 , where /i is an inverse length scale. Thus, Vk/w k ~ /i 2_2fe 
and hence the Vk/w k have a negative naive dimension which decreases drastically with increasing k. This leads to the 
conclusion that the Vk are highly irrelevant couplings. Though irrelevant, one must not set Vk = in calculating the 
noise exponents. We expand the scaling function /; yielding 



Cf ((x,x');t,™>0) 
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with being expansion coefficients depending on x, x', and r. It is important to recognize that = because 
the corresponding terms are not generated in the perturbation calculation. The first term on the right hand side 
of Eq. (^) gives the leading behavior. Thus, vanishes upon setting vi = and we cannot gain any further 
information about . In other words, the vi are dangerously irrelevant in investigating the critical properties of the 
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Our renormalization group improved perturbation calculation comprises the diagrams A to L given in Ref. 
Since they are irrelevant, the couplings proportional to vi have to be treated by inserting [fll 
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into these diagrams. Here </>(p, A ) denotes the Fourier transform of (phi, ). In Schwinger parametrization the 
computation of the diagrams (including insertions) involves summations of the structure 
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We replace these summations by integrations which can be carried out by completing the squares in the exponential. 
In the limit D — > we obtain 

- ({s}f ViKi (!)+••• , (16) 

where the Cj ({s}) are homogeneous functions of the Schwinger paramters Sj of degree zero. After all, the total 
contribution of a diagram can be written as 



I (p 2 , a) = -ViKt (a) jf E S J c i ^ s » 2 ' D (P 2 ' + 
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Here I? (p 2 , {s}) stands for the integrand one obtains upon Schwinger parametrization of the corresponding diagram 
in the usual (f> 3 theory. 

The ellipsis in Eq. ( |l7| ) stands for primitive divergences corresponding to all operators Of' of the generic form 

" 2a 

A p 2 4> n having the same or a lower naive dimension than O"' (a + b + n < I + 2). This myriad of newly generated 
operators is required as counterterms in the Hamiltonian. Inserting either of the O^p , however, does not generate 
C?(". Thus, the renormalization scheme is given by 

o\L =^o^ n + E^ ^L (is) 

i 

^Lc =E45 Een , (19) 

i 

and one solely needs Z^ in calculating the scaling index of Oren ■ Therefore, we refer to o£L as master and to the 
of\ ea as his servants. 

Now we briefly illustrate our perturbation calculation at one-loop level. Upon inserting into the conducting 
propagators of diagrams A and B one finds in e-expansion 



Aoo -2B 0(2 , =g^ T -^i ~—v 2 K 2 (A) +i-v 2 ( A ] + ~ \ p 2 \, (20) 



where G e = (4ir)^ d ^ 2 ^ T (1 + e/2). We learn, that not only primitive divergencies proportional to K 2 ^A^, but also 
" 2 (~ 2 \ 2 

proportional to p 2 A and A I are generated. Thus, one has, at least in principle, to consider O^ 2 ' in liaison with 



0f ] = -\v 2 J d d P (t^j 4> (p, A) (-P, - 



(21) 



where ^ Ai=0 and Pi = 0. The diagrams one obtains upon insertion of these operators are depicted in Fig. [lj. A 
straightforward calculation reveals that neither of the diagrams in Fig. [I] contains primitive divergencies proportional 
to K 2 (^X^j . Hence, these diagrams can be neglected in calculating the scaling index of 0^ 2 \ 



4 




(2) (2) (2) 

FIG. 1. One-loop diagrams obtained by inserting (1) 0\ , (2) 0\ and (3) 0\ respectively. 

The situation is different if one wishes to compute corrections to scaling associated with C'j 2 '* JlJ] . All the diagrams 
in Fig. |l| enter the correction exponent and thus one has to compute and diagonalize a 3 x 3 renormalization matrix. 
This was overlooked by Harris and Lubensky |l3f who erroneously neglected diagrams b, c, d and diagram a with the 

(2) 

0\ insertion . 

By employing dimensional regularization and minimal subtraction we proceed with standard techniques of renormal- 
ized field theory JTTQ - We calculate to two-loop order. From the renormalization group equation and dimensional 
analysis we deduce that the correlation function G scales at criticality as 

G(x,x';A) = |x-x'| 2 - d -"{l + W A' |x-x'|^ + ^(a) |x - x'|^ + • ■ ■ } . (22) 

v and rj are the well known critical exponents for percolation [Q . <f> is the resistance exponent Jl5],[l6) , <j> = 1 + e/42 + 
4e 2 /3087+ (e 3 )- For the noise exponents ipi, I > 2, we obtain here 



ipi = 1 + 



c 2 



7(1 + 0(1 + 20 12348 (l + 3 (l + 2/) 3 
x |313 - 672 7 + Z |3327 - 4032 7 - 8z|4 (-389 + 273 7 ) 

+ I [-2076 + IOO87 + I (-881 + 336 7 )] }} 

- 672 (l + 2 (l + 20 2 *(l + 2/)| + O(e 3 ) . (23) 
7 = 0.5772... denotes Euler's constant and \1/ stands for the Digamma function. With Eq. (pf) the desired scaling 

(n) 

behavior of C R is now readily derived yielding 

C£° ~ Ix-x'^"/" . (24) 

Our result for the noise exponents is in agreement to first order in e with the one-loop calculation by PHL. We 
point out that Eq. (|2^) can be analytically continued to I = 1 and is in conformity with the result for cf> cited above. 
Analytic continuation of ipi to I = and comparison with the available e-expansion results for Db |jl^,^8|,|o| shows 
that -00 = vDb (5) up to order O (e 3 )- Blumenfeld et al. |L9| proved that ipi is a convex monotonically decreasing 
function of I. Note that our result for ipi captures this feature for reasonable values of e. It reduces to unity in the 
limit I — y 00 as one expects from the relation of "0oo to the fractal dimension of the singly connected (red) bonds [^) , 
ipoo = rfrcd^, and Coniglio's proof [^0| of d lc d = l/v. 

In conclusion, we introduced the concept of master and servant operators and showed that it works consistently as 
a tool to describe the multifractal properties of RRN by renormalized field theory. We presented the premier two-loop 
calculation of a family of multifractal exponents. Our result (|2^) is for dimensions near the upper critical dimension 
6 the most accurate analytic estimate of {ipi} that we know of. It fulfills several consistency checks. 

We showed a one-to-one correspondence of the multifractal moments and the master operators. Though a myriad 
of servant operators is involved in the renormalization of the masters O^ 1 ' the scaling behavior of the Zth multifractal 
moment is governed by only. The situation is different for operators which are irrelevant without being mas- 
ters. The scaling behavior of the related quantities is influenced by the whole bunch of operators generated in the 
perturbation calculation. 

The concept of master and servant operators should have many more applications. Indeed, its applicability might be 
a general feature of multifractal systems. Hence, our concept could prove to be a key in understanding multifractality, 
at least from the standpoint of renormalized field theory. 
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